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Abstract 



Using a QCD inspired model developed by our group for particle production, the Interacting Gluon Model 
(IGM), we have made a systematic analysis of all available data on leading particle spectra. These data 
^ , include diffractive collisions and photoproduction at HERA. With a small number of parameters (essentially 

only the non-perturbative gluon-gluon cross section and the fraction of diffractive events) good agreement 
with data is found. We show that the difference between pion and proton leading spectra is due to their 
£^ , different gluon distributions. We predict a universality in the diffractive leading particle spectra in the large 

■ momentum region, which turns out to be independent of the incident energy and of the projectile type. 

OO 

. PACS number(s): 13.85.Qk,11.55.Jy 

a,: 
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jj ■ 1 Introduction 

a, 

In high energy hadron-hadron collisions the momentum spectra of outgoing particles which have the same 
quantum numbers as the incoming particles, also called leading particle (LP) spectra, have been measured 
already some time ago [Q. Recently new data on pion-proton collisions were released by the EHS/NA22 
collaboration Q in which the spectra of both outcoming leading particles, the pion and the proton, were 
simultaneously measured. Very recently data on leading protons produced in eletron-proton reactions at HERA 
with a c.m.s. energy one order of magnitude higher than in the other above mentioned hadronic experiments 
became available |3| . In the case of photoproduction data can be interpreted in terms of the Vector Dominance 
Model Q and can therefore be considered as data on LP production in vector meson-proton collisions. These 
new measurements of LP spectra both in hadron-hadron and in eletron-proton collisions have renewed the 
interest on the subject, specially because the latter are measured at higher energies and therefore the energy 
dependence of the LP spectra can now be determined. 

It is important to have a very good understanding of these spectra for a number of reasons. They are 
the input for calculations of the LP spectra in hadron-nucleus collisions, which are a fundamental tool in the 
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description of atmospheric cascades initiated by cosmic radiation g. There are several new projects in cosmic 
ray physics including the High Resolution Fly's Eye Project, the Telescope Array Project and the Pierre Auger 
Project H for which a precise knowledge of energy flow (LP spectra and inelasticity distributions) in very high 
energy collisions would be very usefull. 

In a very different scenario, namely in high energy heavy ion collisions at RHIC, it is very important to know 
where the outgoing (leading) baryons are located in momentum space. If the stopping is large they will stay 
in the central rapidity region and affect the dynamics there, generating, for example, a baryon rich equation 
of state. Alternatively, if they populate the fragmentation region, the central (and presumably hot and dense) 
region will be dominated by mesonic degrees of freedom. The composition of the dense matter is therefore 
relevant for the study of quark gluon plasma formation . 

In any case, before modelling p — A or A — A collisions one has to understand properly hadron-hadron 
processes. The LP spectra are also interesting for the study of diffractive reactions, which dominate the large 
Xf region. 

Since LP spectra are measured in reactions with low momentum transfer and go up to large xp values, it 
is clear that the processes in question occur in the non-perturbative domain of QCD. One needs then "QCD 
inspired" models and the most popular are string models, like FRITIOF, VENUS or the Quark Gluon String 
Model (QGSM). Calculation of LP spectra involving these models can be found in refs. § and ||. 

In the framework of the QCD parton model of high energy collisions, leading particles originate from the 
emerging fast partons of the collision debris. There is a large rapidity separation between fast partons and sea 
partons. Fast partons interact rarely with the surrounding wee partons. The interaction between the hadron 
projectile and the target is primarily through wee parton clouds. A fast parton or a coherent configuration of 
fast partons may therefore filter through essentially unaltered. Based on these observations and aiming to study 
p — A collisons, the authors of ref. Q proposed a mechanism for LP production in which the LP spectrum is 
given by the convolution of the parton momentum distribution in the projectile hadron with its corresponding 
fragmentation function into a final leading hadron. This independent fragmentation scheme is, however, not 
supported by leading charm production in pion-nucleus scattering. It fails specially in describing the D~ / D + 
asymmetry. A number of models addressed these data and the conclusion was that valence quark recombination 
is needed. Translated to leading pion or proton production this means that what happens is rather a coalescence 
of valence quarks to form the LP and not an independent fragmentation of a quark or diquark to a pion or a 
nuclcon. Another point is that the coherent configuration formed by the valence quarks may go through the 
target but, due to the strong stopping of the gluon clouds, may be significantly decelerated. This correlation 
between central energy deposition due to gluons and leading particle spectra was shown to be essential for the 
undertanding of leading charm production Q . 

In this work we follow the same general ideas of ref. || but with a different implementation. In particular we 
replace independent fragmentation by valence quark recombination and free leading parton flow by deceleration 
due to "gluon stripping". These ideas are incorporated in the model employed by us, the Interacting Gluon 
Model (IGM), which has been used to study energy flow in non-diffractive reactions JII], [l^] and has been 
recently extended to diffractive processes [|13| and also applied to the recent HERA (photoproduction) data 
both on diffractive mass distributions Jl4|] and leading J/vt spectra |lf|. 

We shall study all measured LP spectra including those measured at HERA. We will find and comment 
universal aspects in the energy flow pattern of all these reactions. Universality means, in the context of the 
IGM, that the underlying dynamics is the same both in diffractive and non-diffractive LP production and both 
in hadron-hadron and photon-hadron processes. 



2 The Interacting Gluon Model 



The Interacting Gluon Model (IGM) was introduced some time ago [[IT] and developed by us recently (fj, [l3| 
and proved out to be quite usefull for the study of energy flow. Since the model has been extensively discussed 
in our previous papers we shall present here only the basic ideas and a few formulas leaving the more detailed 
discussion to the appendix. The main aspect of the IGM, shared with minijet models such as HIJING Jlj|, is the 
assumption that hadron-hadron reactions are dominated by multiple and incoherent parton-parton scatterings. 
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Among these, gluon-gluon scatterings are the most important. At very high energies and large scales this is 
a very good approximation. At not very large energies and lower scales one is already moving towards the 
non-perturbative domain and dealing with soft gluons and the incoherence hypothesis might not be valid. 

The soft gluons involved in the collisions studied here are partly pre-existing inside the hadron and partly 
produced by radiation. Pre-existing soft gluons exhibit the same properties as those studied in lattice QCD 
simulations in the strong coupling regime. According to di Giacomo and collaborators Q , the typical correlation 
length of the soft gluon fields is around 0.2 — 0.3 fm. Since this length is still much smaller than the typical 
hadron size, the gluon fields can, in a first approximation, be treated as uncorrelated. As a consequence [ fL8[ 
the number of (soft) gluon-gluon collisions will follow a Poissonian distribution, which was also used in refs. 



|lq , 19, |20|]. In the case of radiated soft gluons, it was recently shown in [|2l| that gluons produced with small 
transverse momenta are independently emitted from the radiating parton, as QCD coherence supresses their 
showering. Consequently, the multiplicity of low pr gluons follows a Poisson distribution, suggesting that the 
collision number follows a Poisson distribution as well. These facts indicate that in the region where perturbative 
results break down, the independent collision approximation may still be a reasonable one. From the practical 
point of view, it was shown in |Tl| that replacing the Poissonian distribution by a broader one does not affect 
the results significantly as long as some mass scale is introduced to cut off the very low x region. 

In the IGM the two colliding hadrons are represented by valence quarks carrying their quantum numbers 
(charges) plus the accompanying clouds of gluons (which represent also the sea qq pairs and therefore should 
be regarded as effective gluons) . In the course of a collision gluonic clouds interact strongly and form a gluonic 
Central Fireball (CF) located in the central region of the reaction. The valence quarks (plus those gluons which 
did not interact) get excited and form Leading Jets (LJ's) (or Beam Jets) which then populate mainly the 
fragmentation regions of the reaction. 

The valence quarks are therefore spectators and the bulk of the reaction and energy deposition occurs 
because of the gluon-gluon collisions, whose number is Poisson distributed with the mean value given by: 

< /TT - ^ y >) = ° 99 W*) G{x , } GW) e (xV „ K 2 mm) j (1) 



dx'dy' a(s) 

where G"s denote the effective number of gluons from the corresponding projectiles (approximated by the 
respective gluonic structure functions) and a gg and a are the gluon-gluon and hadron-hadron cross sections, 
respectively. In the above expression x' and y' are the fractional momenta of two gluons coming from the 
projectile and from the target whereas K min — m /^/s (with mo being the mass of lightest produced state and 
y/s the total c.m.s. energy). Each gluon-gluon collision produces a minifireball (MF). Depending on the energy 
many g — g collisions may happen and energy fractions x = J^. riix\ and y = n,y^ from the target and 
projectile may be deposited in the central region. The probability for depositing the energy fractions x and y 
can be analitically computed in the IGM and is given by the function x(x, y) derived in the appendix. 

Leading particles can also be diffractively produced. In the IGM diffractive dissociation (DD) can be included 
in a simple way, by just requiring that one of the colliding hadrons looses only a very small fraction of its initial 
energy momentum. This was done by imposing cuts on the first moments of the u> function In our earlier 
calculations of energy flow with the IGM we were not concerned with the fragmentation region. The main 
interest was the energy deposition in the central region, which is highly relevant for quark gluon plasma physics. 
Therefore our first LP spectra were in reasonable agreement with data but in the region x^ > 0.8 the model 
prediction was below the experimental points because we had no diffractive component whereras data did not 
discriminate between diffractive and non-diffractive events. Later on we have included diffraction in the IGM 
JlJ. The resulting picture is the same as for non-diffractive events, except that in DD the gluon cloud of the 
projectile interacts only with a subset of the target gluon cloud, which carries a small momentum fraction and 
which we call "Pomeron" (IP). Our Pomeron is essentially just a kinematical restriction which forces one of 
the colliding hadrons to loose only a very small ammount of its energy. Using the diffractive IGM we have 
obtained a good description of the diffractive mass spectra in hadronic diffractive collisions |l3| and also in 
photon-proton collisions |Q. Although it is possible to reformulate the model in the impact parameter space 
and associate diffractive events with peripheral collisions, as done in ref. p2[ |, we prefer here to explore further 
the kinematical interpretation of diffraction. Moreover, in the present formulation we correctly reproduce the 
diffractive peak, which was not obtained in ref. [Q. Apart from the recent EHS/NA22 and HERA LP spectra, 
a new experiment now under consideration at Fermilab will, among other things, address the question of leading 
particles and DD component in the near future p3|. In view of these facts we shall recalculate LP spectra in 
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our model with a DD component properly included and extend it to e — p collisions. 



3 Leading Hadron Spectra 



In the IGM a typical non-diffractive event is shown in Figure la). Colliding particles loose energy fractions x 
and y, forming leading particles with xp = 1 — x and xp = 1 — y. In the figure, V stands for vector meson, 
used later for photoproduction at HERA. We shall consider the reactions p + p —* p + X, tt + + p — > ir + + X, 
K + + p — > K + + X and tt + + p — > ir + + p + X. Later we also address photon-proton reactions in the VDM 
approach: p + V — > p + X. Diffractive processes are illustrated in Figures lb) and lc), where the Pomeron is 
emitted from the target and from the projectile respectively. In Figs, la), lb) and lc) the probability to form 
a CF with mass M = ^/xys is called % (x, y), Xi( x > y) an d xt( x > V) respectively. These functions were derived 
in our earlier works (cf., for example, [|13| ). 

With the functions x( x > y) we obtain the corresponding LP spectra just by changing variables and by adding 
the resulting distributions with proper weights. 

In the lower legs of Fig. 1, leading particles emerge from the collision keeping momentum fraction xl with 
distribution F^pix^) given by: 



F L p(x L ) = dx 

Jo Jo 



dy 



(l-a) X nd (x,y)+ J2 a jX d j{x,v) 

.7=1,2 

(M LP + m Q f 



■ S(x L -l + y)e[xy Si 



= {I- a) dx X nd (x;y = l-x L 



) + ^ otj dxx d j{x;y = 1 - x L ), 

j = l,2 J"min 



(2) 



where a is the total fraction of diffractive events and ou and a2 are the fractions of diffractive events with a 
Pomeron emitted from the upper and lower leg in Fig. 1, respectively. They satisfy the condition a\ + ai = a. 
In the above expression M^p and too denote the mass of the LP and the mass of the lightest CF produced and 
the limits of integration are defined by 

[ ml (M LP + m ) 2 ] 
Xmin = Max T . (3) 

(1 - XL)S s 

The main physical quantities in (0) are the functions (cf. appendix for details of derivation and justification) 



x nd (x,y) = 



with 



Xo 



■>'<j 



■ exp 



2D, 



[(y 2 )(x - {x)f + (x*)(y - {y)f - 2(xy)(x (x))(y (y))} , 



(4) 



and 



x d (x,y) = 



2tt\ D 



D xv = {x 2 )(y 2 } - (xy) 



(x n y m ) 



l r i 

dx'x' n / dy'y' m uj nd {x',y') 
n Jo 



(5) 
(6) 



■ exp 



2D U) 



[( y ])(x - (x 3 )f + (x]){y - ( yj )f - 2(x jVj )(x - (x 3 )){y - (%•))] , (7) 
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with 



D { £ = (*j)(y$) - {x jUj )\ (8) 
{x^yf) = [ Xma *dx'x ,n j V ^dy'y' m u d {x', y '). (9) 



(i 



The index values of j = 1 and j = 2 correspond to the diagrams in Figs, lb) and lc) respectively. Here 
Umax = 1) Umax = y, JmL = x and XmL = 1 and Xo d an d are tne proper normalization factors assuring 
that Jq dxL Flp(xl) = 1. We separately normalize to unity both components entering (|^). This procedure is 
crucial in our case in order to assure the proper overall energy momentum conservation, which is a characteristic 
feature of any implementation of the IGM. 

The weight a is essentially our new parameter. It should be of the order of the ratio between the total 
diffractive and total inelastic (including DD processes) cross sections, &tlt / a tot- It can ' however, differ from 
that ratio due to the different experimental acceptance of DD and non-DD events not considered here. This 
fact makes a a free parameter of the model. We assume it to be independent of the total c.m.s. energy. Indeed 
the fraction of diffractive events with respect to the total number of events or the ratio between diffractive and 
total inelastic cross sections are quantities which depend weakly on the c.m.s. energy of the collision. They may 
depend more strongly on detector coverage and acceptance. The reactions that we discuss here occur at c.m.s. 
energies ranging from 14 up to 100 GeV . This is still a relatively moderate variation in the energy. Since we do 
not address diffraction at the SppS (energies of 200 - 900 GeV) nor at the Tevatron (energy of 1800 GeV) we 
keep our fraction of diffractive events as a constant. This quantity is the only free parameter in this paper. 

The spectral functions uj nd and uj d are the same as in @ and @ [Uj for non-diffractive and for diffractive 



processes, respectively. We refer the reader to |12|, |13|, |14| for details concerning the exact values and character of 
the relevant parameters which were fixed from other previous applications of the IGM. In the case of diffractive 
hadron-proton scattering the function G(y) in eq. (0) represents the momentum distribution of the gluons 
belonging to the proton subset called Pomeron and y is the momentum fraction of the proton carried by one of 
these gluons. We shall therefore use the notation G(y) = Gp(y). This function should not be confused with 
the momentum distribution of the gluons inside the Pomeron, f g /p(P) (see below). 

The moments (q n ), q = x,y (we only require n = 1, 2) are the only places where dynamical quantities like 
the gluonic and hadronic cross sections appear in the IGM. In DD we are selecting a special class of events and 
therefore we must choose the correct dynamical inputs in the present situation, namely Gp{y) and the hadronic 
cross section a appearing in uj d . 

As already mentioned, the Pomeron for us is just a collection of gluons which form a color singlet and belong 
to the diffracted hadron. In early works we have assumed that these gluons behave like all other ordinary gluons 
in the proton and have therefore the same momentum distribution. The only difference is the momentum sum 
rule, which for the gluons in IP is 

l 

dyyG F (y)=p (10) 

where p ~ 0.05 (see |H| and below ) instead of p ~ 0.5, which holds for the entire gluon population in the 
proton. 

In principle, since the Pomeron can not be considered an ordinary particle, one can not define precisely a 
momentum sum rule in the usual sense p4| , p5[ . Nevertheless, in our definition, the Pomeron is really just a 
collection of gluons. As such, when emitted, it carries a momentum fraction, p, of the parent proton. Of course, 
p can fluctuate. In the calculations, it always appears divided by the Pomeron-proton cross section, which is 
another poorly known and fluctuating quantity. Taking all these fluctuations into account would just introduce 
more freedom in the model and make calculations more complicated. We avoid it in this paper as we also avoid 
other sources of fluctuations like impact parameter fluctuations. In ref. |22| impact parameter fluctuations were 
introduced in the non- diffractive IGM and the result was, as expected, a smearing of the original curves. The 
overall effect, however, was not very big. We expect the same to be true for the fluctuations in the Pomeron 
sum rule. 

In juj we have treated the Pomeron structure in more detail and addressed the question of its "hardness" 
or "softness". In order to make contact with the analysis performed by HERA experimental groups we have 
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considered two possible momentum distributions for the gluons inside IP, one hard, fg/p((3) and the other soft 
fg/p(P)- Following the (standard) notation of ref. j2(J, {3 is the momentum fraction of the Pomeron carried by 
the gluons and the superscripts h and s denote "hard" and "soft" respectively. Using a standard choice for the 
Pomeron flux factor, fp/ p {xp), where xp is the fraction of the proton momentum carried by the Pomeron and 
noticing that (3 = the distribution Gp(y) needed in eq. (|lj) is then given by the convolution: 

G h p s (y) = f F/P (x P ) f h g f p {^-) (ii) 

In ref. |Q we have also used Gp(y) = 0.3 ^~ v ' , the same expression already used by us before in ref. p3[ . 
As it was shown, this choice corresponds to an intermediate between "soft" and "hard" Pomeron and it shall 
be used here. 

One of the most interesting conclusions of ref. |l4| was that it is possible to learn something about the 
Pomeron profile studying the diffractive mass spectra. Moreover, our analysis suggests that the "soft" Pomeron 
is in conflict with these data. Only with a very unusual choice of parameters a good agreement could be 
recovered. Considering the large ammount of data already described previously by the IGM, this choice was 
extreme and we concluded therefore that the "soft" Pomeron is disfavoured. The same conclusion was found in 
refs. (26| . The fraction of diffracted nucleon momentum, p, allocated specifically to the IP gluonic cluster and 
the hadronic cross section a are both unknown. However, they always appear in to as a ratio ) of parameters 
and different choices are possible. Just in order to make use of the present knowledge about the Pomeron, we 
shall choose 

a(s) = a Fp (s) = a + b In— , (12) 

so 

where so = 1 GeV 2 and a — 2.6 mb and b = 0.01 mb are parameters fixed from our previous Jl3[ systematic data 
analysis. As it can be seen, a(s) turns out to be a very slowly varying function of a/s assuming values between 
2.6 and 3.0 mb, which is a well accepted value for the Pomeron-proton cross section, and p ~ 0.05 (cf. Q). 

As it will be seen, a good undertanding of the systematics of LP production can be obtained in terms of the 
dynamical inputs contained in eq. ( |l|) . With the exception of a all parameters are fixed. In the next section 
we compare our results given by eq. (|2|) with experimental data. 



4 Results and Discussion 



In Figs. 2a), 2b) and 2c) we present our spectra of leading protons, pions and kaons respectively. The dashed 
lines show the contribution of non-diffractive LP production and the solid lines show the effect of adding a non- 
diffractive component. All parameters were fixed previously and the only one to be fixed was a. For simplicity 
we have neglected the second diagram in Fig. 1, because it gives a curve which is very similar in shape to the 
non-diffractive curve. In contrast, the Pomeron emission by the projectile (Fig. lc) produces the diffractive 
peak. We have then chosen a\ = and ai — a = 0.3 in all collision types. As expected, the inclusion of the 
diffractive component flattens considerably the final LP distribution bringing it to a good agreement with the 
available experimental data [Q . In our model there is some room for changes leading to fits with better quality. 
We could, for example, use a prescription for hadronization (as we did before in [fT2| )) giving a more important 
role to it, as done in ref. ||]. In doing this, however, we loose simplicity and the transparency of the physical 
picture, which are the advantages of the IGM. We prefer to keep simplicity and concentrate on the interpretation 
of our results. In first place it is interesting to observe the good agreement between our curve and data for 
protons (Fig. la) in the low x region. The observed protons could have been also centrally produced, i.e., they 
could come from the CF. However we fit data without the CF contribution. This suggests, as expected, that 
all the protons in this x range are leading, i.e., they come from valence quark recombination. In Figs, lb) and 
lc) we observe an excess at low x. This is so because pions and kaons are light they can more easily be created 
from the sea (centrally produced). Our distributions come only from the Leading Jet (LJ) and consequently 
pass below the data points. A closer look into the three dashed lines in Fig. 2 shows that pion and kaon spectra 
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are softer than the proton one. The former peak at x ~ 0.56 while the latter peaks at x ~ 0.62. In the IGM 
this can be understood as follows. The energy fraction that goes to the central fireball, K = y/xy, is controled 
by the behaviour of the function x(x,y) nd , which is approximately a double gaussian in the variables x and y, 
as it can be seen in expression (0). The quantities (x) and (y) play the role of central values of this gaussian. 
Consequently when (x) or (y) increases, this means that the energy deposition from the upper or lower leg (in 
Fig. 1) increases respectively. The quantities (x) and (y) are the moments of the ui function and are directly 
proportional to the gluon distribution functions in the projectile and target and inversely proportional to the 
target-projectile inelastic cross section. In the calculations, there are two changes when we go from p — p to 
7r — p: 

(i) The first is that we replace cff , by cr^el which is smaller. This leads to an overall increase of the energy 
deposition. There are some indications that this is really the case and the inelasticity in tt — p is larger 
than in p — p collisions . 

(it) The second and most interesting and important change is that we replace one gluon distribution in the 
proton G p (y) by the corresponding distribution in the pion G 7T (y). We know that G p (y) ~ (1 — y) 5 /y 
whereas G n (y) ~ (1 — y) 2 /y, i.e., that gluons in pions are harder than in protons. This introduces an 
asymmetry in the moments (x) and (y), making the latter significantly larger. 

As a consequence, pions will be more stopped and will emerge from the collision with a softer x spectrum. 
This can already be seen in the data points of Fig. 2. However since these points contain particles produced 
by other mechanisms, such as central and diffractive production, it is not yet possible to draw firm conclusions. 

The analysis of the moments (x) and (y) can also be done for the diffractive process shown in Fig. lc). 
Because of the cuts in the integrations in eq. (g), they will depend on xl = 1 — y. We calculate them for 
p+p — > p + X and 7r +p — > 7r + X reactions. For low xl they assume very similar values as in the non-diffractive 
case. For large xl however we find that (x) p ~ (x)^ and (y) p ~ (y)^- The reason for these approximate 
equalities is that in diffractive processes we cut the large y' region and this is precisely where the pion and 
the proton would differ, since only for large y are G P p(y) ~ (1 — y) 5 /y and Gp{y) ~ (1 — y) 2 /y significantly 
different. In rcf. [[u| we have shown that the introduction of the above metioned cuts drastically reduces the 
energy (\/s) dependence of the diffractive mass distributions leading, in particular, to the approximate \ jM\ 
behaviour for all values of yfs from ISR to Tevatron energies. Here these cuts produce another type of scaling, 
which may be called "projectile scaling" or "projectile universality of the diffractive peak" and which means 
that for large enough xl the diffractive peak is the same for all projectiles. The corresponding \ d functions 
will be the same for protons and pions in this region. The cross section appearing in the denominator of the 
moments will, in this case, be the same, i.e., a Fp . The only remaining difference between them, their different 
gluonic distributions, is in this region cut off. This may be regarded as a prediction of the IGM. Experimentally 
this may be difficult to check since one would need a large number of points in large x l region of the leading 
particle spectrum. Data plotted in Fig. 2 neither prove nor disprove this conjecture. The discrepancy observed 
in the proton spectrum is only due to our choice of normalization of the diffractive and non-diffractive curves. 
The peak shapes are similar. 

The EHS/NA22 collaboration provided us with data on ir + + p — > ix + + p + X reactions. In particular they 
present the x distributions of both leading particles, the pion and the proton. Their points for pions and protons 
are shown in Fig. 3a) and 3b) respectively. These points are presumably free from diffractive dissociation. The 
above mentioned asymmetry in pion and proton energy loss emerges clearly, the pions being much slower. The 
proton distribution peaks at ~ 0.6 — 0.8. Our curves (solid lines) reproduce with no free parameter this 
behaviour and we obtain a good agreement with the pion spectrum. Proton data show an excess at large Xf 
that we are not able to reproduce keeping the same values of parameters as before. 

The authors of ref. || tried to fit their measured proton spectrum with the FPJTIOF code and could 
not obtain a good description of data. This indicates that these large x points are a problem for standard 
multiparticle production models as well. In our case, if we change our parameter too from the usual value 
m = 0.35 GeV (solid line) to to = 0.45 GeV (dashed line) we can reproduce most of data points both for 
pions and protons as well. This is not a big change and indicates that the model would be able to accomodate 
this new experimental information. Of course, a definite statement about the subject would require a global 
refitting procedure, which is not our main concern now. 

If, at high energies, the reactions p — p and 7r — p have the same characteristics and if VDM is good 
hypothesis (as it seems to be), then more about the energy flow in meson - p collisions can be learned at HERA. 
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Indeed, as mentioned in |14|, at the HERA electron- proton collider the bulk of the cross section corresponds 
to photoproduction, in which a beam electron is scattered through a very small angle and a quasi-real photon 
interacts with the proton. For such small virtualities the dominant interaction mechanism takes place via 
fluctuation of the photon into a hadronic state which interacts with the proton via the strong force [[| . High 
energy photoproduction exhibits therefore similar characteristics to hadron-hadron interactions. Recent data 
taken by the ZEUS collaboration at HERA Q show that the LP spectra measured in photoproducion and in 
DIS (where Q 2 > AGeV 2 ) are very similar, specially in the large x^ region. This suggests that, as pointed out 
in ]2j|, the QCD hardness scale for particle production in DIS gradually decreases from a (large) Q 2 , which is 
relevant in the photon fragmentation region, to a soft scale in the proton fragmentation region, which is the 
one considered here. We can therefore expect a similarity of the inclusive spectra of the leading protons in high 
energy hadron-proton collisions, discussed above, and in virtual photon-proton collisions. In other words, we 
may say that the photon is neither resolving nor being resolved by the fast emerging protons. This implies 
that these reactions are dominated by some non-perturbative mechanism. This is confirmed by the failure of 
perturbative QCD p0[ , (implemented by the Monte Carlo codes Ariadne and Herwig) when applied to the 
proton frgamentation region. In ref. |29| the LP spectra were studied in the context of meson and Pomeron 
exchanges. Here we use the vector meson dominance hypothesis and describe leading proton production in the 
same way as done for hadron-hadron collisions. The only change is that now we have p — p instead of p — p 
collisions. Whereas this may be generally true for photoproduction, it remains an approximation for DIS, valid 
in the large xl region. 

In Fig. 1 we show schematically the IGM picture of a photon-proton collision. According to it, during the 
interaction the photon is converted into a hadronic (mesonic) state and then interacts with the incoming proton. 
This hadronic state is called V in the upper legs of Figs, la), lb) and lc). At HERA only collisions V — p 
are relevant. The meson-proton interaction follows then the usual IGM picture, namely: the valence quarks 
fly through essentially undisturbed whereas the gluonic clouds of both projectiles interact strongly with each 
other. The state V looses fraction x of its original momentum and gets excited carrying a, xp = I — x fraction 
of the initial momentum. The proton, which we shall call here the diffracted proton, looses only a fraction y of 
its momentum but otherwise remains intact. We shall assume here, for simplicity, that the vector meson is a 
p° and take G p (x) = G n (x) in eq. ([!]). In Fig. 4 we compare our results with ZEUS data. The agreement is 
again good. 



5 Conclusions 

We have analyzed leading particle spectra in terms of the IGM, which includes now also a contribution coming 
from the diffractive processes. The new component improved dramatically the agreement with all existing data 
on hadron-hadron collisions. 

As long as the energy flow is concerned the IGM works extremely well with essentially two parameters: the 
non-perturbative gluon-gluon cross section and the fraction of diffractive events. This should enlarge consider- 
ably its range of applicability in analyses of cosmic ray data. 

At the same time, assuming VDM, we were able to describe equally well the leading proton spectra in e — p 
reactions. Also here the inclusion of a diffractive component provided by the new version of the IGM turns out 
to be crucial to get good agreement with data. 

We have shown that the difference between pion and proton leading spectra is due to their different gluon 
distributions ||3l|| . We predict a universality in the diffractive leading particle spectra in the large momentum 
region, which turns out to be independent of the incident energy and of the projectile type. 



Appendix 

We shall collect together here the main points concerning the Interacting Gluon Model (IGM) scattered through- 
out the literature p0||-|l5[. The IGM is based on the idea that since about half of a hadron momentum is carried 
by gluons and since gluons interact more strongly than quarks, during a collision there is a separation of con- 
stituents. Valence quarks tend to be fast forming leading particles whereas gluons tend to be stopped in the 
central rapidity region. It belongs therefore to the class of models exploring the concept of partons and of hard 
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and semihard collisions (like those presented by Gaisser and Stanev [^9|, Sjostrand p0| , Wang |16| or Geiger 
IH). The latter are collisions between partons at a moderate scale (Q 2 ~ (2GeV) 2 ), which, however, still 
allows for the use of perturbative QCD. The scattered partons form the so-called minijets. At </s = 540 GeV 
the minijet cross section is already 25% of the total inelastic cross section. However, apart from some ambi- 
guity in choosing the semihard scale, these models have to face the problem that even at very high energies 
a significant part of a hadronic collision occurs at scales lower than the semihard one. The atitude taken in 
HIJING jt6), in the Parton Cascade Model |Q and also in the IGM is to extrapolate these quantities to lower 
scales. These extrapolations can be continuously improved, especially in view of the advance of our knowledge 
on non-perturbative effects. There are, for example, models for distribution functions which work at scales as 
low as 0.3 GeV 2 |33|]. As for a one can compute non-perturbative effects in the context of an operator product 
expansion. Inspite of these limitations these models have the advantage of dealing with partons and being thus 
prepared to incorporate perturbative QCD in a natural way. This is welcome since perturbative processes are 
expected to be increasingly important at higher energies. Compared to the other models mentioned above, 
the IGM is simpler because it is designed to study energy flow and makes no attempt to calculate cross sec- 
tions or to follow hadronization in great detail. This simplifies the calculations and avoids time consuming 
numerical simulations. The most important aspect of the IGM, shared with those models, is the assumption of 
multiple parton-parton incoherent scattering which is implicit in the Poissonian distribution of the number of 
parton-parton collisions (which is also used in refs. 0, [l9|, |ll|) used below. 

The IGM is therefore based on the assumed dominance of hadronic collisions by gluonic interactions and 
can be summarized in the following way: 

(i) The two colliding hadrons are represented by valence quarks carrying their quantum numbers (charges) 
plus the accompanying clouds of gluons (which represent also the sea qq pairs and therefore should be 
regarded as effective ones). 

(ii) In the course of a collision the gluonic clouds interact strongly depositing in the central region of the 
reaction fractions x and y of the initial energy-momenta of the respective projectiles in the form of a 
gluonic Central Fireball (CF). 

(Hi) The valence quarks (plus those gluons which did not interact) get excited and form Leading Jets (LJ's) 
which then populate mainly the fragmentation regions of the reaction. 

The fraction of energy stored in the CF is therefore equal to K — ^fxy and its rapidity is Y — | ln^ . These 
two quantities provide then a sort of dynamically calculated initial conditions for any statistical model of 
multiparticle production and that was one of the initial aims of the IGM. 

According to the IGM the CF consists of minifireballs (MF) formed from pairs of colliding gluons. In 
collisions at higher scales a MF is the same as a pair of minijets or jets. In the study of energy flow the details 
of fragmentation and hadron production are not important. Most of the MF's will be in the central region and 
we assume that they coalesce forming the CF. The collisions leading to MF's occur at different energy scales 
given by Q 2 = Xi yi s, where the index i labels a particular kinematic configuration where the gluon from the 
projectile has momentum Xi and the gluon from the target has j/j. We have to choose the scale where we start 
to use perturbative QCD. Many studies in the literature converge to the value Q 2 nin = Pr min — (2.3 GeV) 2 . 
Below this value we have to assume that we can still talk about individual soft gluons and due to the short 
correlation length (found in lattice QCD calculations) between them they still interact mostly pairwise. In this 
region we can no longer use the distribution functions extracted from DIS nor the perturbative elementary cross 
sections. 

The central quantity in the IGM is the probability to form a CF carrying momentum fractions x and y of 
two colliding hadrons. It is defined as the sum over an undefined number n of MF's: 

X(x,y) = ^ 6[x-mxi niXi]8[y - niyi m ■ P(n x ) • • • P(rii) 



\ " 



n ( ai ) 
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The delta functions in the above expression garantee energy momentum conservation and P(rii) is the 
probability to have rij collisions between gluons with xi and yi. If the collisions are independent P(rii) is given 
by: 

(ni) ni exp(-n t ) 



P(ni) 



(A2) 



Inserting Pirii) in (Al) and using the following integral representations for the delta functions: 



E 



fli Xi 



1 f +oc , 

= — / at exp 



fli Xi 



(A3) 



i 

2^ 



du exp 



we can perform all summations and products arriving at: 
\ r+°° 



x{x,y) 



/-t-OQ f-t-OO 
(it / (iw exp[i(te + uy)] exp < j ? 
- oo J — oo 



-i{txi+uyi) 



Taking now the continuum limit: 



we obtain: 



du 



M dx'dy' 
dx' dy' 



(27T) 



(ft 



exp[z(tx + uy)] x 



(A4) 



(A5) 



(A6) 



•1 rl 

x exp / dx 1 i dy'uj(x',y') 



D — i(tx'+uy') 



where 



u(x',y') 



dn 



dx' dy' 



(A7) 



(A8) 



This function lo{x' ', y'), sometimes called the spectral function, represents the average number of gluon-gluon 
collisions as a function of x' e y' , contains all the dynamical inputs of the model and has the form: 



, y = a gg (x'y's) G{xl)G{y > )Q ^ yl _ K 2 ) 
a(s) 



(A9) 



where G"s denote the effective number of gluons from the corresponding projectiles (approximated by the 
respective gluonic structure functions) and a gg and a are the gluon-gluon and hadron-hadron cross sections, 
respectively. In the above expression x' and y' are the fractional momenta of two gluons coming from the 
projectile and from the target whereas K m i n = mo/y/s, with mo being the mass of lightest produced state and 
y/s the total c.m.s. energy, mo is a parameter of the model. 

The integral in the second line of eq. ( A7) is dominated by the low x' and y' region. Considering the singular 
behaviour of the G{x) distributions at the origin we make the following approximation: 



-i(tx' +uy') 



-i(tx' + uy') - - (tx' + uy') 2 



(AlO) 



With this approximation it is possible to perform the integrations in (A7) and obtain the final expression 
for x( x >y) discussed in the main text: 



■ exp 



Xo 



2-KyJ D X y 

1 



2D, 



[(y 2 )(x - (x)f + (x 2 )(y (y)f 2(xy)(x - (x))(y - (y))] 



(All) 
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where 



and 



D xy = (x 2 )(y 2 ) - (xy)'' 



(x n y m ) = j'dxx- [ dyy m oj(x,y), 



Xo is a normalization factor denned by the condition: 



1 ,1 

dx / dyx(x,y)0(xy - K min ) = 1 
o Jo 



(A12) 



(A13) 



In order to evaluate the distribution (All) we need to choose the value of mo, the semihard scale pTmin 
and define G(x) and a gg in both interaction regimes. We take PTmin = 2.3 GeV and mo = 0.35 GeV. These 
are the two scales present in the model. The semihard gluon-gluon cross section is taken, at order a 2 , to be: 



Vqafr, y, s) 



16 



[a s (Q 2 



where 



and 



H = 36 T 



51 AT 3A 2 T 



4xy 8x 2 y 2 



9A 

xy 



In 



H 



1 — T 
T+T 



(A14) 



(A15) 



T = 



i-A 

xy 



X 



4*4, 



(A16) 



The parameter k is the one frequently used to incorporate higher corrections in a s and is 1.1 < K < 2.5 according 
to the choice of G(x), of the scale Q 2 and prmin- For PTmin = 2.3 GeF, k = 2.5. 
The coupling constant is given by: 



a s (Q 2 ) 



12tt 



(33 - 2N f ) In 



(A17) 



where A = 0.2 GeV and Nf = 3 is the number of active flavors. As usual in minijet physics we choose 
Q 2 = p\ min and use the distributions G(x,Q 2 ) parametrized in literature. 

When the invariant energy of the gluon pair s is the interval m\ < s = xys < 4p|, min we are outside 
the pcrturbative domain. Parton-parton cross sections in the non-perturbative regime have been parametrized 
in |}4j leading to a successfull quark-gluon model for elastic and diffractive scattering. Recently these non- 
perturbative cross sections have been calculated in the stochastic vacuum model |35|. The obtained cross 
sections are functions of the gluon condensate and of the gluon field correlation length, both quantities extracted 
from lattice QCD calculations. In order to keep our treatment simple we shall adopt the older parametrization 
for the gluon-gluon cross section used in |M : 



xys 



(A18) 



where cto is the second parameter of the model fjj 
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Figure Captions 



Fig. 1 The schematic description of IGM. In Fig. la) we show a non-diffractive event. The upper (lower) leg 
represents the leading particle with the momentum fraction xp = 1 — x (xl = 1 — y)- In Fig. lb) 
the hadron in the upper leg sends a IP with momentum fraction x and remains a leading particle with 
momentum fraction xp — 1 — x. Fig. lc) corresponds to the case where IP is emitted from hadron in the 
lower leg, which keeps momentum xp = 1 — V- 

Fig. 2 Comparison of our spectra F(xl) with data from ref. jj] for a) leading protons, b) leading pions and c) 
leading kaons. Dashed and solid lines show the non-difractive component and the total curve respectively. 

Fig. 3 a) Comparison of our spectra F{xp) for leading pions with data from ref. [g] in the reaction n + +p — > 
7r + + p + X. Solid and dashed lines correspond to the choices too = 0.35 GeV and mo = 0.45 GeV 
respectively, b) the same as a) for the leading proton spectrum F{xf) measured in the same reaction. 

Fig. 4 Comparison between our calculation and the LP spectrum measured at HERA by the ZEUS collaboration, 
ref. 0. 
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